Predictions for nit and Mw in Minimal Super symmetric Models 

O. Buchmueller^, R. Cavanaugh•^^^ A. Dc Roeck'^''=, J.R. EUis^, H. Flacher^, S. HeinemeyerS, 
G. Isidori^-', K.A. OliveJ, F.J. Ronga'", G. WeigW 

'^High Energy Physics Group, Imperial College, Blackett Lab., Prince Consort Road, 
LondonSW7 2AZ,UK 

^Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, Illinois 60510, USA 
■^Physics Department, University of Illinois at Chicago, Chicago, Illinois 60607-7059, USA 
'iCERN, CH-1211 Geneve 23, Switzerland 
"Antwerp University, B-2610 Wilrijk, Belgium 

^Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA 

sinstituto de Ffsica de Cantabria (CSIC-UC), E-39005 Santander, Spain 

^'INFN, Laboratori Nazionah di Frascati, Via E. Fermi 40, 1-00044 Frascati, Italy 

'Institute for Advanced Study, T. U. Miinchen, Arcisstrafie 21, D-80333 Miinchen, Germany 

■iWilliam I. Fine Theoretical Physics Institute, Univ. of Minnesota, Minneapolis, Minnesota 55455, USA 

•^Institute for Particle Physics, ETH Ziirich, CH-8093 Ziirich, Switzerland 

'DESY, Notkestrasse 85, D-22603 Hamburg, Germany 

Using a frequentist analysis of experimental constraints within two versions of the minimal supersymmetric 
extension of the Standard Model, we derive the predictions for the top quark mass, mt, and the W boson mass, 
Mw- We find that the supersymmetric predictions for both mt and Mw, obtained by incorporating all the 
relevant experimental information and state-of-the-art theoretical predictions, are highly compatible with the 
experimental values with small remaining uncertainties, yielding an improvement compared to the case of the 
Standard Model. 

CERN-PH-TH/2009-223, DESY 09-207, FTPI-MINN-09/44, UMN-TH-2827/09 



One of the most impressive successes of the 
Standard Model (SM) has been the accurate pre- 
diction of the mass of the top quark obtained 
from a fit to precision electroweak measurements 
at LEP and the SLC [T], which agrees very well 
with the value measured at the Tevatron [2]. To 
this may be added the equally successful predic- 
tion of the W mass |l|3j . The successes of these 
comparisons between theory and experiment re- 



quire the incorporation of higher-order quantum 
corrections. In the SM these receive contribu- 
tions from the postulated Higgs boson. Indeed, 
the precision data favour a relatively light Higgs 
boson weighing ^150 GeV [Tj. 

One theoretical framework that predicts such a 
light Higgs boson is supersymmetry (SUSY) [4], 
which also possesses the ability to render more 
natural the electroweak mass hierarchy, contains 
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a plausible candidate for astrophysical dark mat- 
ter, facilitates grand unification, and offers a 
possible explanation of the apparent discrepancy 
between the experimental measurement of the 
anomalous magnetic moment of the muon, {g — 
2)^, and the theoretical value calculated within 
the SM. There have been many analyses of the 
possible masses of particles within the minimal 
supersymmetric extension of the Standard Model 
(MSSM), taking into account the experimental, 
phenomenological and astrophysical constraints. 
For example, we have presented sparticle mass 
predictions [5|6|7j on the basis of a frequentist 
analyses of the relevant constraints in the context 
of simple models for SUSY breaking such as the 
CMSSM (in which the input scalar masses mp, 
gaugino masses mi/2 and soft trilinear parame- 
ters ^0 are each universal at the GUT scale) and 
the NUHMl (in which a common SUSY-brcaking 
contribution to the Higgs masses is allowed to be 
non-universal). For an extensive list of references, 
see [7]. 

These analyses favour relatively light masses 
for the sparticles, indicating significant sensitiv- 
ity of the precision observables to quantum ef- 
fects of supersymmetric particles. It is therefore 
desirable to revisit the successful predictions of 
the SM, in particular the show-case predictions 
of mt and Mw, to see how they are affected in 
the CMSSM and NUHMl. In particular, one may 
ask whether the SM prediction of mt and Mw is 
improved, relaxed or otherwise altered in these 
models. The answer to this key question is non- 
trivial, since low-mass sparticles such as the t and 
b may contribute significantly to the prediction 
of electroweak observables [S], and the (lightest) 
Higgs mass is no longer an independent quantity, 
but also depends on the sparticle masses as we 
discuss below. 

In this Letter, we make supersymmetric pre- 
dictions for mt and Mw within the same frame- 
work as our previous frequentist analyses of the 
CMSSM and NUHMl parameter spaces [STeTZ] . 
The treatments of the experimental, phenomeno- 
logical and astrophysical constraints are nearly 
identical with those in [7]. Here, we employ the 
updated SM value of {g — 2)^ which includes a 
new set of low-energy e+e~ data [9]. The new 



value of {g — 2)^ [10] does not significantly alter 
the regions of the CMSSM and NUHMl parame- 
ter spaces favoured in our previous analyses. 

Our statistical treatment of the CMSSM and 
NUHMl makes use of a large sample of points 
(about 3 X 10^) in the SUSY parameter spaces 
obtained with the Markov Chain Monte Carlo 
(MCMC) technique [IT]. Our analysis is en- 
tirely frequentist, and avoids any ambiguity 
associated with the choices of Bayesian pri- 
ors. The evaluations are performed using the 
MasterCode |5|6|7|12j . which includes the fol- 
lowing theoretical codes. For the RGE running 
of the soft SUSY-breaking parameters, it uses 
SoftSUSY [13], which is combined consistently 
with the codes used for the various low-energy 
observables: FeynHiggs |14|15|16] is used for the 
evaluation of the Higgs masses and a^^^^ (see 
also |17|8j ). for the other electroweak precision 
data we have included a code based on [18ll9j . 
SuFla |20I21| and Superlso |22I23| are used 
for flavour-related observables, and for dark- 
matter-related observables MicrOMEGAs j24j and 
DarkSUSY [25] are used. In the combination of 
the various codes, MasterCode makes extensive 
use of the SUSY Les Houches Accord |26I27| . 

In the SM, the precision of the confrontation 
between theory and experiment is often expressed 
in the {mt,Mw) plane. The experimental values 
of these quantities are essentially uncorrclated [1] 

US], 

mfP = 173.1 ± 1.3 GeV, (1) 
M^^ = 80.399 ± 0.023 GeV, (2) 

shown in Fig. [T] as the black ellipse. In the SM, 
mt is an independent input parameter, whereas 
the relation between the gauge boson masses M\y 
and Mz can be predicted with high precision in 
terms of mt, the Higgs mass, M™, and other 
model parameters, see |28| and references therein. 
The correlation between mt and the prediction 
for Mw is displayed in Fig. [1] (foliated by lines of 
constant Higgs mass, 

A fit of the SM parameters to precision observ- 
ables, e.g., those measured at the Z peak [29] . 
yields indirect predictions for mt and Mf^, and 
hence also a prediction for Mw- The SM predic- 



3 



tion for rrit without including the experimental 
limits on and ex- or including the experi- 

mental measurement of Af^/ is [Tj 

^fit,SM,excLAA, ^ 172.6^1^-3 QcV , (3) 
^fit,SM.incl.Mw ^ 179.3 + 11^6 g^y ^ 

and the SM prediction for Mw, excluding the ex- 
perimental measurement of Mw but either ex- 
or including the experimental measurement of rrit 

ism 

^fit,SM,incl.mt ^ ^ ^ ^^^^ ^ ^.g^ 

The regions of the {mt,My/) plane favoured at 
the 68% C.L. by direct experimental measure- 
ments (HI E]) and in the SM fit jSl E]), shown as 
the dark (blue) contour [3^ in Fig. [T] have signif- 
icant overlap, representing a non-trivial success 
for the SM at the quantum level. However, we 
note that the overlap between the 68% C.L. con- 
tours happens in the region of Higgs mass values 
that are below the exclusion bound from the LEP 
SM Higgs searches, Mf^ > 114.4 GeV [31], in- 
dicating a certain tension between the precision 
observables and the Higgs limit. Indeed, the ex- 
perimental central value of Myy would be reached 
for a Higgs mass as low as ^ 60 GeV. Com- 

bining the indirect measurements, rrit and M\y, 
the best-fit value of M™ ~ 87 GeV, and the 
95% C.L. upper hmit is - 157 GeV [T]. The 

direct searches at the Tevatron currently exclude 
a range 163 GeV < Mf^ < i66 GeV [12], as also 
indicated by a white line in Fig. [U so that the 
range 115 GeV < Mf^ < 150 GeV is favoured 
in a global fit to the SM (including experimental 
bounds) at the 95% C.L. [33]. 

Turning now to our analysis in the case of su- 
persymmetry, we note that the prediction for Myy 
as a function of mt in the unconstrained MSSM 
gives rise to a band in Fig. [1] (shaded green) 
which has only little overlap (shaded blue) with 
the band showing the range of SM predictions for 
Higgs masses above the search limit from LEP. 
This is because the contribution of light super- 
symmetric particles tends to increase the pre- 
dicted value of Mw compared to the SM case. 
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Figure 1. The 68% C.L. regions in the {mt,Mw) 
plane predicted by a SM fit excluding the LEP 
Higgs constraint fj/, and by CMSSM and NUHMl 
fits including the LEP Higgs mass constraint, 
compared with the experimental measurements 
from LEP2 and the Tevatron shown as the black 
ellipse. The medium gray (red) and the dark 
(blue ) shaded regions show the SM prediction, fo- 
liated by lines of constant values. The light 
gray (green) and the dark (blue) regions show the 
prediction of the unconstrained MSSM rang- 
ing from light to heavy SUSY particles. 



Furthermore, the overlap region (corresponding 
to the situation where all supersymmetric par- 
ticles are heavy) is limited because, in contrast 
to the SM, the value of Mh is not an indepen- 
dent parameter in the MSSM, but is calculable in 
terms of the sparticle masses with an upper limit 
~ 135 GeV [IS]. 

We have performed fits in the CMSSM and the 
NUHMl including all relevant experimental in- 
formation as specified in [7], i.e., we include all 
precision observables used in the SM fit shown in 
Fig. [1] (except T^/, which has a minor impact) as 
well as constraints from (g — 2)^, flavour physics, 
the cold dark matter (CDM) relic density and 
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the direct searches for the Higgs boson and su- 
persymmetric particles. The direct experimental 
measurements of A/^y and , on the other hand, 
have not been included in these global fits. The 
results of our fits in the CMSSM and the NUHMl 
are also displayed as 68% C.L. contours in Fig.[Tl 
and show remarkably good agreement with the 
experimental measurements of Mw and mt- 

The 68 and 95% C.L. regions in the {nit, M\y) 
plane found in the CMSSM (NUHMl) fit are 
shown in more detail in the upper (lower) panel 
of Fig. [2 The fits within the MSSM differ from 
the SM fit in various ways. First, the number 
of free parameters is substantially larger in the 
MSSM, even restricting ourselves to the CMSSM 
and the NUHMl. On the other hand, more ob- 
servables are included in the fits, providing ex- 
tra constraints. We recall that in the SM fits 
{g — 2)^ and the i3-physics observables have a 
minor impact on the best-fit regions, and are not 
included in the results shown above, which are 
taken from [1] (see e.g. [34] for an alternative 
approach), while the relic density of CDM can- 
not be accommodated in the SM. Furthermore, 
as already noted, whereas the light Higgs boson 
mass is a free parameter in the SM, it is a func- 
tion of the other parameters in the CMSSM and 
NUHMl. In this way, for example, the masses 
of the scalar tops and bottoms enter not only di- 
rectly into the prediction of the various observ- 
ables, but also indirectly via their impact on M/j. 
This provides additional motivation for including 
the experimental constraints on into the fits 
in the MSSM. 

In Fig. [HI we show the results of the same fit as 
in Fig. [21 but now in the {Mji,mt) plane for the 
CMSSM (NUHMl) in the upper (lower) panel. 
The LEP lower limit of 114 GeV is applicable 
in the CMSSM |35l36j . but cannot always be di- 
rectly applied in the NUHMl, since there are re- 
gions of the NUHMl parameter space where the 
hZZ coupling is suppressed relative to its value in 
the SM [37]. We use the prescription given in [7] 
to calculate the x^(M/i) contribution for points 
with suppressed hZZ couplings, and see in the 
lower panel of Fig. [3] a significant set of NUHMl 
points with Mh <C 114 GeV: these reflect the 
shape of the Ax^ function in the right panel of 
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Figure 2. The 68% and 95% C.L. regions in 
the {rat, M^w) planes for the CMSSM (upper) and 
for the NUHMl (lower), for fits that do not in- 
clude the direct measurements of mt and Mw , but 
do incorporate the appropriate LEP constraint on 



Fig. 4 of [7] . The appearance of a local minimum 
at Mh ~ 100 GeV in the lower plot of Fig. [3] is 
statistically not significant. It sensitively depends 
on the details of the implementation of the Higgs 
search bounds. 

We now turn to the single- variable functions 
for mt and Mw- In the upper panel of Fig.[4l we 
show the functions for mt in the CMSSM and 
NUHMl as solid and dashed lines respectively 
with Mw included in the fit (as before, the di- 
rect measurement of mt is not included in this 
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Figure 3. The 68% and 95% C.L. regions in the 
{Mh,mt) planes for the CMSSM (upper plot) and 
for the NUHMl (lower plot), for fits that do not 
include the direct measurements of mt and Mw , 
but do incorporate the appropriate LEP constraint 
on Mfi- 



fit). Comparing the results with the SM fit, we 
find that these rise more sharply, in particular for 
larger values of m^, than they would in the SM fit, 
indicating that the upper bound on mt from the 
indirect prediction in the MSSM is significantly 
reduced compared to the SM case. We find the 
68% C.L. ranges 



fit. CMSSM, incl.A4"w 

fit.NUHMl.incl. 

rUt ' 



173. 8^i^ ? GcV, 
169.5l^ j GcV. 



(7) 

(8) 



lower central values for mt in both the CMSSM 
and NUHMl in better agreement with the exper- 
imental result ([T]). The reduction in the upper 
bound on mt reflects, in particular, the fact that 
the additional contribution from the t and b en- 
ters with the same sign as the leading SM-type 
contribution to the precision observables that is 
proportional to TOj. A non- vanishing contribu- 
tion from superpartners therefore tends to reduce 
the preferred value of mt compared to the SM fit. 
It should be noted in this context that the smaller 
uncertainties in mt found in the supersymmetric 
fits compared to the SM case (particularly in the 
CMSSM) can in part also be attributed to the fact 
that a larger set of observables has been used in 
the CMSSM and NUHMl fits. 

For the W boson mass, we find the func- 
tions including mt in the fit in the CMSSM (solid) 
and NUHMl (dahed) shown in the lower panel of 
Fig. 31 and the corresponding 68% C.L. ranges 



Af, 
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W 

fit, NUHMl, incl.mt 



W 



= 80.379j:°:gji^ GeV, 
= 80.370t°:Jj?^ GeV. 



(9) 
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Comparing with the SM fit result ([3]), we find 



The best-fit values of these predictions are sub- 
stantially higher than the SM prediction © based 
on precision electroweak data (in particular in the 
CMSSM) and are closer to the experimental value 
([2]). again with smaller uncertainties. 

We summarize our main results in Fig. [5l The 
upper (lower) panel compares the experimental 
measurement of mt (Mw) with the predictions of 
a SM fit to precision electroweak data and our fi- 
nal predictions in the CMSSM and NUHMl. The 
resulting agreement of the final predictions for mt 
with the experimental value (Qp is remarkable, al- 
most embarrassingly good in the CMSSM case, 
and very good in the NUHMl. Compared to the 
SM fit, the best-fit values for Mw in the CMSSM 
and NUHMl are closer to the experimental value 
([2]), and in the CMSSM case the best-fit value 
lies within the experimental 68% C.L. range. We 
conclude that the CMSSM and NUHMl pass with 
flying colours the test of reproducing the success- 
ful SM predictions of mt and Mw, even improv- 
ing on them. We can only hope that this probe 
of SUSY at the loop level will soon be made even 
more precise with the discovery of sparticles at 
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Figure 4. The functions for rrit (upper panel) 
in the CMSSM (solid) and NUHMl (dashed) ex- 
cluding the direct mt mass measurement but in- 
cluding all the other experimental information. 
The corresponding functions for Mw (lower 
panel) excluding the direct Mw mass measure- 
ment but again including all the other experimen- 
tal information. 
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Figure 5. The 68% C.L. ranges for mt (upper 
panel) and Mw (lower panel) including (from top 
to bottom) the experimental average, and the pre- 
dictions of the SM ( not incl. the Alf^ limits ) f]]/, 
CMSSM and NUHMl fits, using all the available 
information except the direct mass measurement. 



